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Transformations 

Executive  Summary 


This  program  had  two  objectives  that  were  pursued  simultaneously.  These 
objectives  were:  (i)  development  of  process  models  designed  to  produce  a  new  family  of  tough 
and  novel  ceramic  materials  with  specific  microstructures  and  (iii)  coordinate  the  finding  of  the 
proposed  research  with  the  ongoing  parallel  efforts  on  the  material  characterization  addressed 
from  the  atomistic  and  continuum  length  scales,  where  correlations  are  made  between  atomic 
structure  and  bonding.  This  final  report  summarizes  the  findings  of  PI  Dr.  A.  Sayir  of  Case 
Western  Reserve  University  and  contains  new  information  relating  development  of  process 
models  designed  to  produce  a  new  mullite  single  crystal.  It  was  concluded  that  the  growth  of 
single  crystal  mullite  is  controlled  by  kinetics  and  it  is  a  non-equilibrium  process.  Liquid-liquid 
immiscibility  and  and  slow  diffusion  rate  of  cations  did  not  permit  the  growth  of  3:2  mullite 
single  crystals.  We  investigated  through  numerical  modeling  effort  for  solidification  of 
representative  important  oxide  crystals,  2Al203*Si02  mullite,  by  the  vertical  Bridgman  technique. 
The  interactive  effects  of  radiative  heat  transfer  and  convection  on  solidification  of  2Al203»Si02 
was  investigated  experimentally  and  isolated  numerically  in  order  to  optimize  the  growth 
conditions  for  solidification  of  the  materials.  A  combined  conduction-convection-radiation 
model  of  solidification  was  developed  for  the  oxide  crystals,  which  are  transparent  to  radiation 
below  wavelength  of  6  //m,  and  opaque  to  radiation  in  the  rest  of  the  spectrum. 

Mullite  exhibits  remarkable  thermo-mecahnical  stability  Crystal  grown  from  melt  show 
a  higher  degree  of  oxygen  vacancy  ordering  and  resistance  to  dislocation  motion.  Columns  of  Oc 
are  not  randomly  distributed  along  [3-10].  Vacancy  rich  colums  tend  to  form  ordered  domains 
with  super  lattice  dimensions  (planes  ±1/3,  ±1,  ]A  and  ±2/3, ±2,  !4).  Vacancy  rich  colums  are 
frequently  seen  along  [131]  and  [261],  Simplistic  model  proposed  and  compared  with  the 
electron  optic  results.  Aluminum  and  silicon  occupied  tetrahedral  site  whereas  A1  also  had 
octahedral  site  occupancy.  Model  did  not  reflect  the  loss  in  oxygen.  It  was  concluded  that  the 
modeling  of  the  structure  of  2Al203*Si02  mullite  required  refined  simulations  using  supercell  to 
describe  the  periodicity  of  vacancies. 

In  a  parallel  study  Co-PI  Prof.  W.  Kriven  of  University  of  Illinois  at  Urbana-Champaign 
investigated  martensite  toughnening  mechanisms.  Toughening  mechanisms  include 
transformation  toughening  (e.g.,  with  zirconia,  AV  =  positive),  and  transformation  weakening 
causing  debonding  of  interphases  (e.g.,  cristobalite,  enstatite,  AV  =  negative).  The  role  of  unit 
cell  shape  changes  remains  relatively  unexplored,  and  can  lead  to  ferroelastic  domain 
rearrangements,  large  force  actuation  or  shape  memory  ceramics.  Basic  science  research  was 
conducted  to  identify,  and  crystallographically  characterize,  potentially  new  martensitic  phase 
transformations  in  yttrium  niobate  (YNbCL)  and  dysprosium  titanate  (DyTiOs),  as  examples  of 
model  systems  for  academic  understanding,  as  well  as  potentially  technologically  useful 
ceramics.  A  unique,  newly-built,  quadrupole  lamp  furnace  capable  of  temperature  over  2400  °C 
in  air,  in  conjunction  with  synchrotron  radiation  enabled  in  situ  studies. 


Nomenclature 


A  =  surface  area,  (m2) 

Bi  =  Biot  number,  (hR  /  KJ 

C  =  specific  heat,  (J  /  kg  K) 

F  =  radiation  function 

Gr  =  Grashof  number,  (g0BR>Th  /  v2) 

K  =  thermal  conductivity,  (W  /  m  K) 

=  unit  norma!  vector 

N%  -  radiation-conduction  number,  (oTh3R  /  Ks) 

N ,  =  crucible  radiation-conduction  number,  (oTh3R  /  Ka) 


Pr 

q 

R 

T 

* 


=  Prandtl  number,  (v/a) 

=  dimensionless  radiative  flux,  {q  foT*  ) 
=  radius,  (m) 

=  temperature,  (K) 


=  Velocity  Vector, 


<  V 


(Gr)2 


Greek 

X,  Xc  =  wavelength,  cut-off  wavelength,  (microns) 
p  =  dynamic  viscosity,  (N  sec  /  m2) 
v  =  Kinematic  viscosity,  (N  sec  /  m2) 

p  =  density,  (kg  /  m3) 

o  =  Stefan  Boltzmann  constant,  (W  /  m2  K4) 

0  =  dimensionless  Temperature  (T  /  Th) 

P  =  expansion  coefficient,  (1  /  K) 

e  =  emissivity 


Subscripts 

a  =  crucible 

c  =  cold  zone 

h  =  hot  zone 

I  =  melting  point 

i  =  crystal 

r  =  radial 

z  =  axial 

°o  =  ambient  parameter 

Radiation  View  Factors 
WW  =  wall-to-wall  radiation  exchange  factor 
WS  =  wall-to-interface  radiation  exchange  factor 
SW  =  interface-to-wall  radiation  exchange  factor 
SS  =  interface-to-interface  radiation  exchange  factor 


1.  INTRODUCTION 

The  brittleness  of  ceramics  for 
structural  applications  continues  to  be  a 
difficult  and  unsolved  problem, 
particularly  at  high  temperatures  and  in 
oxidizing  environments.  The  high 
specific  strength  and  stiffness  of 
ceramic  materials  offers  great  promise 
in  high  temperature  applications. 

Research  to  improve  reliability,  either 
through  composite  or  bulk  toughening, 
has  yet  to  develop  viable  competitors  to 
the  high  temperature  metallic  super¬ 
alloys.  However,  an  increase  in 
fracture  toughness  could  make  ceramic 
matrix  composites  a  desirable 
alternative.  There  are  no  alternative  materials  to  ceramics  for  high  temperature  applications.  At 
present,  the  interface  between  engineering  design  and  materials  science  currently  consists  of  a 
collection  of  empirically  developed  and  marginally  verified  material  properties  obtained  by  ad 
hoc  experiments.  These  are  laborious,  time  consuming  and  expensive.  They  could  however,  be 
accelerated  and  much  for  efficiently  carried  out  with  the  help  of  information  gained  from  model 
experiments  providing  information  about  the  relevant  features  of  each  single  component.  In 
particular,  if  one  understand  how  a  phase  transformation  occurs  and  under  what  conditions  it  can 
be  manipulated,  a  ceramic  materials  designer  would  be  an  a  much  better  position  to  expediently 
design  a  composite  which  takes  advantage  of  the  phase  transformation.  We  selected 
2Alj03*Si02  mullite  as  a  case  study.  The  motivation  for  the  selection  of  mullite  was  related  to 
its  exceptional  deformation  resistance  at  very  high 
temperatures.  Mullite  has  a  very  complex 
structure  as  shown  in  Fig.  1.  Yet,  the  exceptional 
resistance  to  plastic  deformation  at  temperatures 
exceeding  1500  °C  can  not  be  explained  by  the 
complexity  of  the  structure.  The  resistance  to 
deformation  does  not  scale  with  the  melting 
tmperature  as  is  the  case  for  most  oxides.  The 
anticipated  theoretical  defct  chemistry  reactions 
are  shown  in  Fig.  2  for  x  values  0,25  and  0.4, 
corresponding  to  3:2  and  2:1  mullite. 

The  performance  of  this  class  of  materials  depends  directly  on  the  quality  of  the  crystals 
used,  which,  to  a  large  extent,  is  determined  by  the  uniformity  of  temperature  and  concentration 
at  the  growth  interface  during  processing.  Consequently,  heat  and  mass  transfer  to  and  from  the 
interface  during  growth,  and  their  relationship  to  the  thermal  environment  established  by  the 
furnace  and  the  fluid  flow  in  the  melt,  play  key  roles  in  controlling  the  solidification  process. 
There  has  been  a  tremendous  drive  in  the  past  decade  to  optimize  and  improve  the  conditions, 
configurations,  and  procedures  for  growing  and  processing  crystals  but  there  has  been  only  a  few 
experimental  and  modeling  work  for  this  family  of  materials.  We  investigated  through 


Mullite:  2Si4+  -*2AI3*  +  V"0 

AI2(AI2+2xSi2_2X)Oio_x 

3AI203-2Si02  2A!202-Si02 

x=0.25  x=0.4 

Fig.  2  Chemistry  of  mullite 


Temperature  field  Concentration  field  Velocity  field 

Fig.  3  Fusion  of  material  modeling  and  computational  tools.  Temperature,  concentration 
and  velocity  fields  for  marangoni  number  2.4  and  Raynolds  number  24. 


numerical  modeling  effort  the  solidification  of  2Al203»Si02  mullite,  by  the  vertical  Bridgman 
technique. 


Status  of  Process  Model:  Combined  Heat  Transfer  and  Fluid  Flow  Analysis  of  Semi- 

Transparent  Crystal  Solidification  - 

The  processing  of  the  2Al203»Si02  mullite  required  temperatures  of  2000  °C  and  greater 
and  heat  transfer  at  these  very  high  temperatures  are  dominated  by  the  radiation  heat  transfer 
mode.  Yet,  most  of  the  numerical  models  developed  for  solidification  of  crystals  have  neglected 
the  effects  of  radiation  heat  transfer  within  the  crucible.  Brown  (1988)  has  made  a 
comprehensive  review  of  the  papers.  In  general,  these  studies  apply  to  solidification  of 
semiconductors  and  cannot  be  generalized  to  represent  the  behavior  of  oxide  crystals.  The 
earliest  investigations  of  the  effect  of  radiation  on  nonopaque  crystals  is  due  to  Viskanta  and  co¬ 
workers  (Tarshis  et  al.;  1969  and  O'Hara  et  ah;  1968,  Viskanta;  1975).  In  these  studies,  the 
formidable  general  problem  is  reduced  to  a  one-dimensional  analysis  of  steady  state  conduction 
and  radiation  heat  transfer  as  it  occurs  between  two  partially  transparent  media  in  intimate 
thermal  contact.  The  models  implemented  in  these  works  neglect  many  of  the  complexities  of 
the  actual  growth  process.  Nevertheless,  they  show  that  the  temperature  distributions  and  the  net 
thermal  fluxes  in  a  growing  semi-transparent  crystal  are  drastically  modified  by  the  participation 
of  the  media  in  the  radiation  transfer  process.  These  investigations  also  conclude  that  radiative 
heat  transfer  stabilizes  the  interface.  This  occurs  because  the  temperature  gradients  at  the 
interface  are  steeper  in  combined  radiative-conductive  systems  than  for  the  purely  conductive 
cases.  This  conclusion  is  later  disputed  by  Abrams  and  Viskanta  (1974).  This  finding  is 
corroborated  by  two  papers  due  to  Antonov  et  al.(  1 980)  and  Yuferev  and  Vasil'ev  (1987).  These 
two  investigations  present  a  more  realistic  treatment  of  conduction  and  radiation  heat  transfer  for 
growth  of  sapphire  by  the  Stepanov's  method.  The  analyses  are  still  one-dimensional  and  still 
the  radiation  heat  transfer  can  lead  to  instability  in  the  crystallization  front.  It  is  also  postulated 
that  the  pronounced  stretching  of  the  solid-liquid  interface  observed  experimentally  may  be  due 
to  radiative  transfer.  Interaction  of  multi-dimensional  radiation  exchange  with  convection  and 


conduction  was  studied  by  Kassemi  et  al  (1989,  1990)  in  the  context  of  crystal  growth  by  vapor 
transport.  They  used  a  finite  difference  control  volume  approach  to  discretize  the  diffusive  and 
convective  terms  and  an  element-to-node  zonal  approach  to  represent  the  radiative  exchange. 

The  results  show  that  through  interaction  with  convection  radiation  significantly  affects 
the  crystal  shape.  It  was  also  shown  that  surface  and  internal  radiation  exchange  have  opposing 
effects  on  the  interface  curvature.  The  first  studies  to  focus  on  solidification  of  semi-transparent 
oxides  from  melt  in  realistic  crystal  growth  configuration  have  been  due  to  Brandon  and  Derby 
(1991,  1992).  Lin  and  Motakef  (1993)  considered  solidification  of  BSO  in  a  vertical  Bridgman 
configuration.  They  solved  the  combined  radiation  conduction  convection  problem  indirectly  by 
iterating  between  the  finite  element  code  ABAQUS  to  solve  the  radiation-conduction  problem  in 
the  solid  and  the  fluids  code  FIDAP  to  solve  the  conduction-convection  problem  in  the  melt.  In 
summary,  the  interactive  effects  of  radiative  heat  transfer  and  convection  on  solidification  of 
YNb04  and  2Al203#Si02  will  be  investigated  experimentally  and  isolated  numerically  in  order  to 
optimize  the  growth  conditions  for  solidification  of  both  materials. 


DESIGN  OF  PROCESS  MODEL 
Physical  Model:  Mathematical  Formulation 

Oxide  crystals  have  been  usually  grown  by  the  Czochralski  technique.  In  the  Czochralski 
process,  heat  is  supplied  radially  from  the  side  to  the  crucible  containing  the  charge  and  removed 
axially  at  the  solidification  front.  This  sets  up  an  unsteady  three-dimensional  destabilizing 
buoyancy-driven  convection  that  produces  periodic  fluctuations  of  the  interface  and  generation 
of  rotational  impurities.  Since  the  requirement  for  uniformity  is  more  stringent  for  oxide  crystals 
(with  length  scales  on  the  order  of  centimeters)  than  for  semiconductors,  researchers  have  been 
forced  to  exploit  other  techniques  for  growing  these  materials.  The  vertical  Bridgman  and  float 
zone  techniques  have  been  a  favorite  alternative  and  has  been  used  for  the  investigation  of 
2Al203»Si02  mullite  solidification  to  reveal  the  intrinsic  capabilities  of  these  materials.  In 
vertical  Bridgman  technique  the  crucible  is  aligned  with  the  gravitational  vector,  convection  in 
the  melt  is  quite  weak  and  the  axisymmetric  heat  input  and  extraction  from  the  sides  results  in  a 
more  stable  axisymmetric  solid-liquid  interface. 


A  cross-sectional  view  of  the  cylindrical  ampoule  in  a  typical  vertical  Bridgman  furnace 
configuration  the  ampoule  which  encapsulates  both  the  solid  and  the  melt  is  pulled  downwards  in 
the  gravitational  field  through  the  temperature  profile  maintained  along  the  bore  of  the  furnace, 
at  a  constant  velocity,  Upull.  The  three-zone  furnace  consisted  of  a  hot  section  separated  from  the 
cold  section  by  an  insulated  region.  The  temperatures  in  the  hot  and  cold  zones  were  maintained 
by  the  insulating  walls,  therefore,  they  werea  assumed  to  be  uniform  at  temperatures  Th  and  Tc, 
respectively.  The  insulated  region  was  not  actively  controlled;  thus,  the  temperature  of  the  bore 
surface  in  this  region  varies  linearly  between  the  hot  and  cold  ends.  As  a  result  the  bore  surface 
temperature  can  be  assumed  by  the  piecewise  continuous  profile.  The  boundary  conditions  has 
to  be  studied  accurately  for  the  2A!203»Si02  mullite  crystals.  The  details  of  the  geometry, 
temperature  levels,  and  thermophysical  and  radiative  properties  are  different  for  each  crystal. 
Thus,  the  details  of  the  thermophysical  and  radiative  properties  were  investigated  experimentally 


and  isolated  numerically  in  order  to  optimize  the  growth  conditions  for  solidification  of  both 
2Al203»Si02  mullite. 


Radiation  heat  transfer 

In  solidification  of  semiconductors  usually  both  the  crystal  and  the  melt  are  opaque  to 
thermal  radiation  and  both  phases  have  relatively  high  thermal  conductivities.  Therefore, 
conduction  is  the  dominant  mode  of  heat  transfer  in  both  the  solid  and  the  melt.  Oxide  crystals, 
on  the  other  hand,  are  usually  semi-transparent  to  thermal  radiation  in  the  solid  phase  and  almost 
opaque  in  the  melt.  They  also  have  relatively  low  thermal  conductivities  in  both  phases. 
Therefore,  heat  transfer  during  the  solidification  of  oxide  crystal  is  governed  by  an  intricate 
balance  between  convection  and  conduction  in  the  melt  and  conduction  and  radiation  in  the 
solid.  In  a  sense  the  solid  acts  basically  like  a  light  pipe  through  which  the  interface  looses 
considerable  amount  of  heat  (by  emission)  to  the  cold  sections  of  the  crucible  wall  or  directly  to 
the  furnace,  if  the  crucible  is  also  transparent.  Therefore,  radiation  heat  transfer  plays  dominant 
role  in  shaping  the  crystal  interface.  Through  its  interaction  with  natural  convection,  it  can  also 
drastically  modify  the  flow  fields  in  the  melt  and  alter  the  segregation  patterns  at  the  growth 
front. 

Governing  Equations 

The  heat  transfer  mechanisms  which  govern  solidification  of  oxides  are  complicated. 
There  was  a  strong  heat  transfer  link  through  radiation  and  natural  convection  between  the 
furnace  and  the  crucible.  Consequently,  the  longitudinal  temperature  distribution  of  the  bore 
surface  was  approximately  imposed  on  the  outer  surface  of  the  crucible.  Heat  transfer  inside  the 
crucible  was  dominated  by  conduction  and  convection  in  the  opaque  melt  and  by  conduction  and 
radiation  through  the  semitransparent  crystal.  Therefore,  at  the  inner  crucible  surface,  in  the  melt 
region,  there  is  a  balance  between  the  conductive  fluxes.  Similarly,  at  the  inner  crucible  surface, 
in  the  crystal,  there  is  a  balance  among  the  radial  conductive  fluxes  through  the  ampoule  and  the 
crystal  and  the  net  radiative  flux  at  the  wall  surface.  Since  the  pulling  velocity  of  the  ampoule 
was  small  a  quasi-steady  approach  was  adopted  and  the  negligible  convective  terms  which  arise 
due  to  the  ampoule  motion  were  not  considered.  The  final  location  and  shape  of  the  interface  (at 
quasi-steady-state)  is  of  course  determined  by  a  balance  among  the  conductive  fluxes  on  both 
sides  of  the  solidification  front,  the  net  radiative  flux  at  the  interface,  and  the  released  heat  of 
fusion.  For  the  materials  considered  here,  and  again  due  to  the  very  slow  pulling  speed,  the 
latent  heat  of  fusion  is  negligible  compared  to  the  other  terms  in  the  energy  balance  but  it  is  still 
retained  in  the  interface  equation  for  sake  of  completeness.  Furthermore,  since  it  was  assumed 
that  the  longitudinal  axis  of  the  ampoule  is  aligned  with  the  gravitational  vector,  the  bore 
temperature  profiles  established  by  the  furnace  are  circumferentially  uniform  and,  therefore, 
axisymmetric  conditions  needs  to  be  exploited.  Following  the  quasi-steady  approach,  the  energy 
equation  describing  heat  transfer  in  the  melt,  the  crystal  and  the  ampoule  is  written  as  follows. 

1 .  Crystal 

oo 

V20  +  {Nrs  )f(v  •  qri  =  0  (EQ  1) 

0 


2.  Melt 


3.  Ampoule 
V20  =  O 

These  equations  are  subject  to  the  following  thermal  boundary  conditions: 


(EQ3) 


2.  Melt-Ampoule  Interface: 
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3.  Crystal-Ampoule  Interface: 
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3.  Ampoule-Air-Furnace  Interface: 
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(EQ4) 


(EQ5) 


(EQ  6) 


4.  Ampoule  Top  and  Bottom: 

9=6C,  Qh  (EQ7) 


In  addition  to  these  boundary  conditions,  the  energy  equation  is  also  subject  to  the  following  in¬ 
terface  conditions. 


0=0„ 

and 


\KS, 


( V0) •  nis  =  -V0 •  n,s-{Nrs)f[qh  •  -  (Sr)(Pe) 


(EQ  8) 


(EQ  9) 


An  example  of  the  proposed  work  can  be  examined  using  commercially  available  crystals  of 
YAG  and  BSO  and  their  radiative  properties.  The  total  divergence  of  radiative  flux  reduces  to 

oo 

j(v  •  qri  jrfA  =  0  (EQ  10) 

0 


The  net  radiative  flux  at  the  wall  is  given  by 


j(^-L)rfA  =  ^F0_A(Ac,  TX 

-feJl-xM"  TX(r')WW(r',  r)dA 

A, 
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Similarly,  the  net  radiative  flux  at  the  interface  is  given  by 


/(^•^)rfA  =  £jF0_At(Ac,  TX 


-fat U(Ae,  TX(r')WS{r',  r)dA 
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(HQ  11) 


(EQ  12) 


The  challenge  for  the  2Al203«Si02  mullite  investigation  was  the  quantitative 
determination  of  Equations  1  through  10  in  a  form  suitable  for  inclusion  2Al203»Si02. 


The  fluid  flow  in  the  melt  which  is  primarily  driven  by  the  buoyancy  force  is  described  by 
the  continuity  and  balance  of  momentum  equation  and  is  subject  to  non-slip  boundary  conditions 
at  all  solid  boundaries. 

1 .  Continuity: 

V«$=0  (EQ  13) 


2.  Momentum: 

(Gr',2)$-Vtf=  -VP  +  V2$(GrU2)(e-8m)  (EQ  14) 


3.  Non-Slip  Boundary  Condition: 

$=0  (EQ  15) 


Finally,  these  equations  are  subject  to  symmetry  conditions  along  the  centerline  of  the  ampoule: 


$  *nr  =  0 

0 


(EQ  16) 


In  writing  the  above  equations  its  was  also  inherently  assumed  that  the  thermophysical 
properties  of  the  solid  and  the  melt  although  different  are  constant  throughout  each  phase,  the 
gap  between  the  furnace  bore  and  the  crucible  wall  is  small  so  that  the  external  radiation 
exchange  between  elements  on  these  two  surfaces  is  basically  one  to  one  (view  factor  is  one), 
and  all  surfaces  in  the  model  are  treated  as  diffuse  for  the  radiation  exchange  calculations.  For 


the  solidification  processes  which  involve  doping,  the  species  transport  equations  and  their 
associated  boundary  conditions  which  include  solutal  rejection  at  the  growth  interface  will  be 
formulated  and  incorporated  in  the  numerical  model.  Solution  of  the  final  system  of  integro- 
partial  differential  equations  will  describe  the  dynamic  coupling  and  temporal  link  between  the 
species,  fluid  flow  temperature  and  radiation  fields. 


Numerical  Formulation  and  Results 


The  finite  element  code  FIDAP  was  adopted  to  solve  the  solidification  problem.  FIDAP 
has  to  be  modified  significantly  in  order  to  incorporate  a  radiative  heat  transfer  model  based  on 
an  exchange  factor  scheme,  similar  to  the  DEF  method  described  by  Naraghi  and  Kassemi 
(1989)  into  the  code.  The  diffusive  and  convective  terms  in  the  governing  differential  equations 
were  discretized  based  on  the  regular  Galerkin  formulation.  Discretization  of  the  integral  terms 
describing  radiation  exchange  using  the  exchange  factor  method  were  necessary  to  achieve  a  set 
of  algebraic  equations  which  can  be  incorporated  in  both  finite  element  and  finite  difference 
codes.  The  exchange  factor  methods  have  been  already  used  with  both  finite  element  and  finite 
difference  techniques  and  have  produced  very  accurate  results  (Kassemi  and  Naraghi  1993, 
Saltiel  and  Naraghi,  1990).  Nine-node  quadratic  continuum  elements  and  three-node  quadratic 
surface  elements  was  used  throughout  the  mesh.  The  radiation  flux  terms  can  then  be  evaluated 
at  these  integration  points. 


The  resulting  set  of  discretized  algebraic  equations  will  be  solved  using  a  segregated 
approach  which  tackles  every  single  degree  of  freedom  individually  (Engleman,  1993).  That  is, 
the  solution  matrix  is  decomposed  into  smaller  loosely  connected  sub-matrices  each  associated 
with  only  one  conservation  equation.  The  smaller  matrices  are  then  solved  in  a  sequential 
manner  using  gaussian  elimination.  Because  these  matrices  are  loosely  connected  the  method 
has  a  slower  convergence  rate  compared  to  the  coupled  direct  solvers.  However,  the  scheme  is 
fairly  robust  because  it  has  a  large  radius  of  convergence  and  does  not  require  a  good  starting 
guess.  The  degrees  of  freedom  are  decoupled,  it  is  much  less  demanding  with  regard  to 
computer  storage  needs  as  compared  to  coupled  solvers. 

The  solidification  of  mullite  was  treated  by  tracking  the  interface  as  a  free  surface.  In 
contrast  to  the  coupled  solvers  where  the  new  interface  position  is  part  of  the  solution  vector,  in 
the  segregated  approach  an  iterative  scheme  is  adopted.  First,  the  conservation  equations  are 
solved  assuming  a  fixed  interface  location  and  using  the  continuity  of  temperature  at  the 
interface  (Eq.  1 1)  as  a  boundary  condition  to  the  energy  equation.  Next,  the  required  change  in 
the  interface  location  is  determined  by  satisfying  the  interface  energy  balance  condition  (Eq.  12). 
The  updated  interface  position  is  then  used  again  to  solve  the  conservation  equations.  This 
iterative  process  is  continued  until  a  converged  solution  is  attained.  During  each  simulation,  as 
the  solidification  front  moves  and  changes  shape,  the  interface  position  is  tracked  and  the 
position  of  the  nodes  are  changed  along  straight  lines  (spines)  in  proportion  to  the  movement  of 
the  interfacial  nodes. 


High  resolution  digital  image 


Fig.  4  Liquid-liquid  immiscibility  of  molten  mullite  liquid  (left).  Strategic  experimentation 
and  interoperatibility  with  existing  sotware  FIDAP  results  shwing  mild  change  of  Si/Al  at  the 
interior  of  the  melt  (right) 


The  results  for  the  temperature-  concentration-  and  velocity  field  are  shown  in  Fig.  3  for 
dimensionless  Marangoni  number  2.4  and  dimensionless  Raynolds  number  24.  Figure  3  is 
results  for  the  large  diameter  solidification  condition.  The  fundamental  challenge  the  coupling 
between  temperature  and  velocity  field  to  produce  concentration  field.  As  will  be  discussed 
further  in  following  section  single  crystal  growth  of  2Al203*Si02  mullite  was  preceded  with 
liquid-liquid  immiscibility.  The  high-resolution  digital  images  of  laser  heated  float  zone  system 
revealed  that  the  melt  separated  into  two  immiscible  liquid  as  determined  from  the  different 
optical  transmission  characteristic,  Fig.  4.  Finite  element  code  FIDAP  results  also  included  in 
Fig.  4  and  reveals  mild  change  of  Si/Al  ratio  in  the  molten  liquid.  The  molten  zone  of  the 
growth  experiments  with  the  3Al203»2Si02  mullite  compositions  were  quenched  to  analyze 
global  chemical  partitioning  in  the  melt.  The  experimental  results  could  not  identify  in 
significant  compositional  differences  the  length  scales  of  molten  zone  within  the  range  40  to  600 
H m.  Finite  element  code  FIDAP  results  showed  spatial  change  in  the  Si/Al  ratio.  Although  this 
change  was  very  small,  it  was  suffient  to  have  different  degree  of  polymerization  of  tetrahedral 
connectivity  of  the  molten  Al203»Si02  composition.  This  had  tremendous  effect  on  the  viscosity 
of  the  melt  and  by  far  most  important  parameter  was  the  temperature  of  the  melt  to  alter  the 


degree  of  the  polymerization.  In  fact,  to  fabricate  single  crystal  2Al203*Si02  mullite,  melt  had  to 
be  heated  at  temperatures  exceeding  2100  °C.  This  means  at  least  200  °C  superheating  to 
produce  single  crystals.  The  starting  composition  3Al203*2Si02  mullite  could  be  melted  at 
temperature  1 800  °C  and  wide  range  microstructures  could  be  produced  at  temperature  region 
between  1800  to  2200  °C. 


Growth  Experiments 


The  crystal  growthinvestigation  was  initiated  to  determine  in  more  detail  the  growth  behavior  of 
mullite  single  crystal  fibers.  The  primary  objective  of  this  work  was  to  determine  the  growth 
characteristics  of  single  crystal  mullite  fibers  grown  by  the  Laser  Heated  Floating  Zone-Method 
(LHFZ)  and  compare  the  results  with  fibers  grown  by  Edge-Defined  Film  Fed  Growth  (EFG). 
The  effect  of  different  supercooling  conditions  and  starting  compositions  were  compared.  Also 
included  in  this  study  are  the  effects  of  melt  modulation  to  minimize  void  formation.  Another 
objective  was  to  study  the  tensile  failure  of  mullite  fibers.  The  results  and  possible  mechanisms 
are  discussed  by  examining  the  relationships  between  growth  characteristics,  fracture  stress,  and 
mullite  fiber  microstructures.  The  growth  characteristic  of  fibers  in  the  SiCVA^Oa  system  has 
not  been  documented  and  findings  presented  here  provide  a  valuable  supplement  to 
understanding  the  tensile  strength  characteristics  of  mullite  filaments. 


Oriented  solidification 


Bubbles  + 

Liquid-Liquid  immiscibility 


Fig.  5  The  effect  of  melt  modulation  on  controlling  void  formation  is  shown  in  a). 
Bubbles  are  restricted  to  regions  below  the  laser  scan  and  do  not  rise  to  the  growth 
interface.  The  phase  separated  region  remains  invariant  with  pulling  direction  upward  a) 
and  downward  b). 


There  are  at  least  two  major  sources  of  heterogeneities  in  aluminosilicate  melts  during  laser 
heated  float  zone  growth.  The  first  heterogeneity  is  the  development  of  gas  bubbles  and  the 


second  is  liquid-liquid  immiscibility.  These  heterogeneties  have  profound  effects  on  the  fiber 
growth  process  and  resulting  mechanical  and  microstructural  fiber  characteristics.  The  effect  of 
bubbles  in  the  high  viscosity  Si02-Al203  melt  on  fiber  growth  can  be  minimized  by  a  melt 
modulation  technique,  a  new  heating  technique  which  eliminates  bubbles  from  the  growth  front 
region.  The  SiC>2-Al203  melt  has  high  a  viscosity  and  is  known  to  entrap  bubbles  coming  from 
the  source  materials.  The  bubbles  can  be  entrapped  in  the  fiber  as  voids  (flaws)  and  can  either  be 
present  on  the  surface  or  in  the  bulk  of  the  material.  The  dynamic  nature  of  the  molten  zone  and 
bubble  formation  or  removal  was  monitored  in  situ  from  high-resolution  digital  images.  Melt 
modulation  effectively  eliminated  void  formation  at  or  near  the  surface  of  solidifying  fibers. 
Figure  6  shows  a  photograph  of  a  mullite  fiber  grown  with  the  melt  modulation  technique  in 
which  the  transverse  scan-length  was  larger  than  the  molten  zone  diameter.  This  technique 
agitates  the  molten  liquid  of  Si02-Al203.  There  was  no  indication  of  bubble  formation  above  the 
laser  scanning  plane  during  growth.  The  absence  of  voids  was  confirmed  by  subsequent  SEM 
and  TEM  examination.  Fibers,  with  very  few  exceptions,  fractured  from  surface  flaws.  The 
growth  of  void  free  fibers  at  high  pull-rates  was  not  possible  with  stationary  heating.  The 
dominant  feature  of  stationary  laser  heating  is  that  a  large  amount  of  energy  is  absorbed  in  a 
small  area  in  a  short  time  and,  hence,  does  not  eliminate  the  void  formation  as  effectively  as  the 
scanning  technique.  The  melt-modulation  technique  has  at  least  two  additional  advantages. 

First,  a  more  uniform  temperature  distribution  can  be  achieved  in  the  molten  zone.  Uniform 
growth  temperature  is  a  prerequisite  for  any  crystal  growth  system  and  it  is  even  more  important 
for  the  growth  of  mullite  single  crystals  which  show  a  strong  tendency  for  facet  formation  (Fig. 
7).  Second,  mixing  in  the  melt  can  be  improved  by  vibration  of  the  melt  with  different  scan 
frequencies. 


Fig.  6  Complex  coupling  between  different  transport  mechanism  at  the  interface.  High 
resolution  digital  images  showing  faceted  growth. 

The  second  major  source  of  heterogeneities  in  SiCVA^C^  melts,  liquid-liquid  immiscibility  in 
the  high  AI2O3  compositions,  has  been  suspected,  but  so  far  has  not  been  unambiguously 
demonstrated.  The  review  of  stored  images  revealed  that  liquid  adjacent  to  single  crystal  mullite 
crystals  was  more  transparent  than  in  the  interior  of  the  melt  (Figs.  1 , 2,  and  3).  The  position  of 
the  darker  dome  shaped  interior  liquid  was  not  effected  by  the  laser  impingement  position. 
Reflected  optical  microscopy  of  longitudinal  cross-sections  of  quenched  molten  zone  specimens 
showed  scattering  intensity  features  that  correspond  with  the  images  shown  in  Figs.l  and  2. 


Hence,  we  interpreted  these  features  arising  from 
an  immiscibility  in  the  liquid,  possibly  due  to 
compositional  and/or  structural  differences  in  the 
melt. 

Usually  systems  with  a  miscibility  gap  in 
the  liquid  phase  have  density  differences  between 
the  two  liquid  phases  large  enough  for  the  two 
phases  to  remain  separated  as  long  as  they  are  kept 
above  the  monotectic  temperature, TM,  and  below 
the  convolution  temperature  Tc.  However,  the 
anticipated  density  differences  in  the  Si02-Al203 
system  are  not  expected  to  be  large  enough  to  form 
separate  layers.  The  tendency  for  layer  separation 
was  examined  by  comparing  upward  and 
downward  growth  conditions.  The  spatial  position 
of  the  darker  area  was  invariant  with  respect  to 
growth  direction  (Fig.6)  and  remained  attached  to 
the  source  rod.  The  density  differences  in  the  two 
immiscible  regions,  therefore,  are  insignificant. 
Yet,  there  was  apparently  enough  difference  in 
surface  tension  between  the  two  liquid  phases  to 
form  separate  regions  in  the  molten  zone.  Similar 
immiscibility  regions  were  observed  for  all 
compositions  studied. 


Fig.  7  The  melt  is  kept  in  its  initial 
undercooled  condition,  a  concave 
solid/liquid  interface  formed  (left). 
Resulting  polycrystalline  mullite,  showing 
paralellel  grains  of  around  10  pm  in  width 
running  over  the  section  (right) 


The  immiscibility  of  Si02-Al203  in  this 
composition  range  was  first  suggested  by  Aramaki 

and  Royin  a  study  in  which  they  stated  that,  on  the  basis  of  Dietzefs  field  strength29  calculations, 
a  wide  (60  mole  %  A1203 ;  3:2  mullite)  miscibility  gap  might  be  expected  at  temperatures  below 
2000  °C.  Galakhov  and  Konolova,  and  McDowell  and  Beall  also  suggested  possible  liquid 
phase  separation.  Recently  Rishbud's,  Rishbud  and  Pask's  and  Kuo's  calculations  using  mullite 


Fig.  8  Square  shaped  grains  separated  by  low  angle  grain  boundaries  (left).  (001)  perpendicular  to 
fiber  axis  (middle).  Elongated  mullite  grains  separated  by  an  alumino-silicate  amorphous  phase. 
Short  range  orering  of  oxygen  vacancies  (right). 


and  alumina  as  chemical  species  in  the  solution  model  predicted  stable  liquid  immiscibility  at 
temperatures  up  to  and  greater  than  2100  °C.  However,  a  miscibility  gap  in  the  Si02-Al203 
system  at  high  AI2O3  concentrations  and  at  T  >  TM  has  not  been  experimentally  observed.  The 
liquid-liquid  miscibility  gap  observations  presented  here  provide  the  first  experimental  evidence 
for  liquid  immiscibility  for  these  compositions  and  temperatures.  Our  experimental  observations 
suggest  that  liquid  immiscibility  remains  a  dominant  heterogeneity  in  the  molten  zone  during 
crystal  growth  at  steady-state  conditions  even  for  higher  temperatures. 

Fibers  grew  as  single  crystal  or  polycrystalline  mullite  as  a  function  of  the  growth 
conditions  and  seeding.  The  polycrystalline  fibers  can  be  grown  from  all  compositions  studied 
in  the  present  work  at  steady-state  conditions  as  long  as  the  melt  stays  in  its  initial  undercooled 
condition,  Fig.  7.  The  SEM  and  TEM  analysis  revealed  that  polycrystalline  fibers  consist  of 
highly  elongated  and  faceted  mullite  crystals  bounded  by  aluminosilicate  glass  (Figs.  8).  Fiber 
cross  sections  are  quite  round  and  faceting  along  the  fiber  length  is  absent.  The  absence  of 
alumina  in  these  fibers  was  confirmed  by  X-RAY  and  TEM  characterization.  These  results  and 
the  microstructure  shown  in  Figs.  8  are  similar  to  those  observed  by  Aksay  and  Pask  and  Kriven 
and  Pask20  on  rapidly  cooled  mullite.  Their  interpretation  of  the  microstructures  on  the  basis  of 
the  supercooling  of  the  liquid  is  in  agreement  with  our  observations. 

The  mullite  fibers  studied  are  transparent,  their  cross  sections  are  not  circular  but  facetted  with 
dimension  of  100  to  150  pim.  They  can  enclose  very  fine  silica  rich  intergranular  phases  as 
shown  by  electron  microprobe  analysis  (EPMA-WDX)  separating  elongated  grains.  These  grains 
have  an  Al/Si  atomic  ratio  of  -4.4  which  corresponds  to  x~0.45  (Al/Si  =  (4+2x)/(2-2x).  The 
growth  direction  corresponds  to  [001], 


Fig.  9  WDX  profile 


The  liquid-liquid  immiscibility  and  its  consequences  on  single  crystal  growth  can  be 
deduced  by  the  observation  of  seeding  and  subsequent  single  crystal  growth  of  mullite.  Fig.7 


Composition 

a,  [nm] 

b,  [nm] 

c,  [nm] 

M70.4(2) 

0.7586 

0.7684 

0.2885 

M71.0(2) 

0.7589 

0.7681 

0.2888 

M72.0(2) 

0.7589 

0.7690 

0.2899 

Table  1  The  lattice  parameters  for  M70.4(2),  M71(2)  and  M72(2)  materials  were  calculated 
from  x-ray  diffractometry  scans. 

Is  a  representative  of  single  crystal  mullite  growth  using  polycrystalline  seed  material.  As 
growth  initiates,  Fig.  6,  the  fibers  are  polycrystalline  and  the  liquid/solid  interface  is 
macroscopically  very  rough  (Fig.  3a).  By  increasing  the  melt  temperature,  a  planar/liquid-solid 
interface,  can  be  established  as  shown  in  Fig.  6.  Thus  single  crystal  growth  can  been  achieved 
by  establishing  steady-state  growth  conditions  from  a  polycrystalline  fiber  interface  (Fig.  5  and 
7).  Although  the  exact  values  are  yet  not  known  for  TM  and  Tc,  the  planar  growth  interface  in 
the  present  work  was  attained  below  Tc,  for  all  compositions  studied.  Thus  single  crystal  growth 
was  achieved  without  a  complete  monotectic  invariant  transition  in  the  liquid  state.  This 
suggests  that  the  faster  diffusion  rates  obtained  by  increasing  temperature  to  achieve  planar 
growth  overcome  the  constitutional  supercooling. 

The  lattice  parameter  differences  measured  over  the  narrow  composition  range  followed 
the  reported  trend  of  a  and  c  increasing  with  alumina  content  (Table  1).  However  the  b 
parameter  did  not  consistently  decrease  as  has  been  reported.  The  differences  were  not 
statistically  significant.  Both  lattice  parameter  values  and  ordering  values  fall  within  the 
reported  range  for  material  compositions  more  alumina  rich  than  2:1  but  clearly  lower  than  5:2. 
As  observed  in  TEM,  mullite  single  crystals  are  free  of  dislocations,  low  angle  boundaries  and 
voids.  Single  crystal  mullite  showed  a  high  degree  of  oxygen  vacancy  ordering  (x  =  0.47,  both  e 
and  f  reflections  present,  Cameron).  Single  crystal  growth  can  also  be  achieved  by  changing  the 
modulation  frequency  (scan-velocity)  of  the  laser  heating.  The  degree  of  supercooling 
diminishes  as  the  solution  agitation  is  increased.  These  observations  lead  us  to  suggest  that 
formation  of  mullite  in  aluminosilicate  melts  is,  in  fact,  preceded  by  immiscibility.  Thus,  the 
mullite  growth  is  accompanied  by  formation  of  a  stagnant  SiC>2-rich  liquid  layer  at  the 
solidifying  interface. 

Figure  6  also  demonstrate  facet  formation  in  single  crystal  mullite  fibers.  The 
crystallographically  induced  metastability  manifests  itself,  as  a  preferential  growth  of  specific 
planes.  This  continued  until  sufficient  departure  from  the  equilibrium  concentration  developed 
and  was  followed  by  sudden  transversing  the  solidifying  interface  direction.  This  formed  a 
"bamboo"  type  fiber  as  shown  in  Figs.  7.  This  was  accompanied  by  a  change  between  the 
equilibrium  contact  angle,  between  the  facet  face  and  the  rest  of  the  interface  (Fig.7)  such  that 
fiber  growth  was  not  centrosymmetric  with  respect  to  the  fiber  axis  and  molten  zone.  The 
selective  adsorption  in  the  solution  causes  a  decrease  of  the  surface  free  energy  of  fast  growing 
plane  and  the  wetting  angle.  The  deviation  of  the  wetting  angle  with  respect  to  average  surface 
energy  of  the  fiber  liquid  system  is  depicted  clearly  in  Figs.  6.  The  growth  of  high  strength 
fibers  requires  tight  control  of  "surface  tension  balance  "  at  the  point  of  intersection  of  the 
interface.  Bolling  and  Tiller  showed  that  this  depends  on  the  temperature  gradient,  surface 
tension,  and  the  entropy  of  fusion  per  unit  volume  for  a  pure  material.  This  implies  that  in  a 


heterogeneous  liquid  phase  in  which  the  activity  of  each  component  varies  with  concentration, 
the  faceting  tendency  will  be  further  complicated  by  mixing  free  energy.  This  suggests  the 
necessity  to  study  the  orientation  relationship  between  supercooling  and  preferred  growth 
direction  of  mullite  single  crystals.  SEM  analysis  revealed  that  the  fibers  were  indeed  strongly 
faceted  (Fig.  6)  and  typically  exhibit  two  large  flat  surfaces  (presumably  {110})  about  the 
circumference.  The  cross-sectional  shape  can  become  highly  oblate  as  the  facet  size  increases. 
Faceting  also  occurs  along  the  fiber  length  causing  sharp  changes  in  fiber  diameter  (bamboo  type 
fiber). 


Mullite  structure 

Mullite  is  a  solid  solution  of  Al203and  Si02  the  chemical  formula  of  which  can  be  written  as 
AlVI2(AlIV2+2xSilv2.  2x)O,0_x.  The  mullite  structure  is  obtained  from  the  silimanite  one  (x=0  , 

Al203  Si02)  by  the  substitution  of  Si4+  by  Al3*  cations.  This  substitution  is  counterbalanced  by  the 
creation  of  oxygen  vacancies  2Si4+  +  O2  ->  2A13+  +  Q.  The  molar  ratio  3Al203-2Si02  or  x=0.25 
corresponds  to  the  stable  form  of  mullite  (also  called  “stoichiometric  mullite”).  However  the  low 
interdiffusion  rates  of  aluminum  and  silicon  through  the  mullite  crystal  make  the  synthesis  of 
metastable  compounds  A1VI2(A1  2+2xSi  2-  2x)O|0.x  possible,  with  a  large  range  of  x.  Composition 
from  around  x=0.25  (3Al203-2Si02)  up  to  x=0.8  have  been  reported.  The  synthesis  route  (i.e. 
solidified  from  the  melt,  obtained  by  solid  state  sintering  from  ceramic  powders  or  sol  precursors 
mixed  at  the  atomic  or  nanometer  scales),  the  rate  of  cooling,  the  presence  or  absence  of  a-Al203 
seeds,  and  the  degree  of  silica  evaporation  are  some  of  the  experimental  factors  explaining  the 
differences  in  the  compositions  and  structures  reported  in  the  literature.  They  have  in  turn  a  direct 
influence  on  the  structural  and  functional  properties  of  mullites.  As  an  example,  the  ordering 
processes  of  the  oxygen  vacancies  have  been  pointed  out  as  a  key  factor  to  understand  the 
remarkable  thermo-mechanical  stability  of  mullite. 


The  structure  of  mullite  is  derived  from  that  of  silimanite  (x=0  ,  Al203  Si02)  in  which  A106 
octahedral  chains  sharing  edges  run  parallel  to  the  c-axis  and  are  cross-linked  by  doublets  of 
Si04  and  A104  tetrahedra.  In  the  silimanite,  A1  and  Si  atoms  are  completely  ordered  in  the 
tetrahedral  chains  with  a  unit  length  a~7.6A,  b~7.7A,  c~5.8A.  The  mullite  structure  is  obtained 
from  the  silimanite  one  by  the  substitution  of  Si4+  by  AI3+  cations.  This  substitution  leads  to  an 
apparent  long  range  disorder  of  A1  and  Si  cations  in  the  tetrahedral  sites.  As  a  result  the  c- 
dimension  for  the  average  structure  becomes  half  that  of  the  silimanite.  The  cation  substitution  in 
the  tetrahedra  is  counterbalanced  by  the  creation  of  oxygen  vacancies  2Si4+  +  O2'  -*  2A13+  +  []. 
Moreover  the  two  cations  which  have  lost  bridging  oxygen  atoms  are  coordinate  to  only  three 
oxygen  atoms,  and  therefore  are  displaced  from  the  T  site  to  a  new  position  IVA1*.  The  adjacent  Oc 
site  is  displaced  off  the  symmetry  center  toward  the  IVA1*  site  to  a  new  Oc*  site  and  is  now  three 
coordinated  by  tetrahedral  cations. 

With  the  assumption  that  the  vacancies  are  randomly  distributed  in  the  structure,  mullite  is 
described  by  an  average  structure  presented  in  Table  I. 


Fig.  10  Atomic  displacements  around  oxygen  vacancy  in  mullite,  average  structure. 
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The  stable  form  of  mullite,  also  called  “stoichiometric”  mullite  corresponds  to  x=0.25  (or 
3Al203-2Si02). 


Oxygen  vacancy  ordering  in  alumina  rich  mullites 

Mullites  show  in  the  reciprocal  space  a  complex  scheme  of  sharp  and  diffuse  satellite  reflections 
and  diffuse  streaks,  in  addition  to  the  Bragg  reflections  expected  from  the  Pbam  space  group, 
Fig.  10.  They  reveal  a  non  random  distribution  of  the  oxygen  vacancy.  These  additional 
reflections  change  their  sharpness  and  position  according  to  the  alumina  content  and  preparation 
mode.  Structural  models  have  been  proposed  in  the  literature  to  explain  the  various  diffraction 
patterns  obtained,  they  correspond  either  to  long  range  ordering  or  short  range  ordering  of  the 
vacancies. 


Sharp  satellite  reflections  are  associated  to  a  long  range  order  of  the  vacancies.  There  are  visible 
for  SAD  pattern  in  the  [010]  zone  axis.  Alumina  rich  mullites  (xs0.4,  2:1  molar  ratio)  show  in 
the  reciprocal  space  Various  structural  models  are  proposed  in  the  literature  to  explain  these 


diffraction  patterns.  Nakajima  et  al1.  report  a  high  degree  of  order  in  the  oxygen  vacancy 
arrangement.  The  diffraction  patterns  and  the  HREM  images  of  alumina  rich  mullites  are 
induced  by  an  incommensurate  modulation  of  the  structure  with  the  formation  of  ordered 
domains  exhibiting  antiphase  relations  parallel  to  the  c-axis,  and  an  incommensurate  period 
along  the  a-axis.  In  opposition  to  this  description,  Rahman  propose  that  these  diffraction  patterns 
are  caused  by  a  short-range  ordering  of  the  oxygen  vacancies.  The  structure  is  then  described  via 
inter-vacancy  correlation  vectors,  the  sequence  and  frequency  of  which  depend  on  the  alumina 
content,  Fig.  11. 

Investigation  of  images  taken  around  Scherzer  defocus  with  a  zone  axis  [310]  in  a  rather  thick 
zone  showed  that  the  corresponding  FFT  had  extra  spots  at  ±h/3,  ±k,  ±1/2  and  diffuse  streaks 
around  these  spots  parallel  to  [001]*,  as  in  a  SAD  pattern.  The  HR  image  shows  variations  in 
intensity  of  the  dots,  Rhaman  have  attributed  these  variations  to  variations  in  vacancy 
concentration  along  the  incident  beam.  Their  image  simulations  for  similar  conditions  of  defocus 
show  that  highest  concentrations  give  brightest  spots.  The  binary  image  has  been  derived  by  a 
thresholding  on  the  gray  level,  to  select  the  brightest  dots  (last  class  of  the  histogram,  2800- 
3200).  These  dots  are  distributed  preferentially  along  specific  directions  as  [131]  and  [261],  rows 
pointed  by  the  arrows.  However  this  directions  of  alignment  do  not  contribute  to  the  contrast 
enhancement  around  ±h/3,  ±k,  ±1/2,  seen  on  the  FFT  of  threshold  image,  and  already  described 
for  the  SAD  patterns.  In  order  to  evidence  the  origin  of  this  periodicity  on  the  vacancy 
distribution,  an  inverse  FFT  images  built  from  respectively  the  spots  a  and  a’,  has  been 
superimposed  to  the  corresponding  threshold  image.  The  resulting  image  reveals  domains  of  a 
few  nanometers  in  which  the  vacancies  tend  to  form  to  form  a  superlattice,  Fig.  12. 

From  this  first  analysis  it  can  be  concluded  that  disparities  in  the  vacancy  distribution  along 
[310]  in  the  Oc  columns  are  seen.  The  vacancy  rich  columns  tend  in  turn  to  form  ordered 
domains  with  super  lattice  dimensions  (planes  ±1/3,  ±1,1/2  and  ±2/3,  ±2,1/2).  In  addition  rows  of 
several  vacancies  rich  columns  are  frequently  seen  along  [131]  an  [261].  The  domains  of  oxygen 
vacancy  ordering  is  shown  in  Fig.  12. 


Fig.  1 1  Extra  reflections  has  been  due  to  a  short  range  ordering  of  oxygen  vacancies  which  are 
characterized  in  terms  of  distribution  frequency  of  intervacancy  vectors. 


1  Y.  Nakajima  et  al.  American  Mineralogist  66  (1981)  142-147 


Fig.  12  Oxygen  vacancy  ordering  in  mullite  as  revealed  first  time. 


CONCLUSIONS 


Directionally  solidified  fibers  with  nominal  mullite  compositions  of  SA^C^SiC^  were  grown 
by  the  laser  heated  float  zone  (LHFZ)  method.  High  resolution  digital  images  from  an  optical 
microscope  evidence  the  formation  of  a  liquid-liquid  miscibility  gap  during  crystal  growth. 
Experimental  evidence  shows  that  the  formation  of  mullite  in  aluminosilicate  melts  is  in  fact 
preceded  by  liquid  immiscibility.  SEM  analysis  revealed  that  the  fibers  were  strongly  faceted 
and  that  the  facets  act  as  critical  flaws.  Examined  in  TEM,  these  mullite  single  crystals  are  free 
of  dislocations,  low  angle  boundaries  and  voids.  Single  crystal  mullite  showed  a  high  degree  of 
oxygen  vacancy  ordering.  Regardless  of  the  starting  composition,  the  degree  of  order  observed 
in  polycrystalline  fibers  was  lower  than  that  observed  in  the  mullite  single  crystals.  Mullite 
formation  (stable  or  metastable  crystalline  phases,  polycrystalline  or  single  crystal)  is  a  function 
of  the  melt  temperature.  Formation  of  a  stagnant  liquid  layer  and  phase  separation  is  strongly 
dependent  on  the  diffusional  properties  of  the  liquid  and  directly  affects  the  type  of  crystal 
grown. 

It  was  proposed  that  the  growth  of  single  crystal  mullite  is  controlled  by  kinetics  and  thus  it  is  a 
non-equilibrium  process.  Liquid-liquid  immiscibility  and  slow  diffusion  rate  of  cations  did  not 
permit  the  growth  of  3:2  mullite  single  crystals.  Mullite  exhibits  remarkable  thermomechanical 
stability.  Crystals  grown  from  the  melt  show  a  higher  degree  of  oxygen  vacancy  ordering  and 
resisance  to  dislocation  motion.  Columns  of  Octend  to  forma  ordered  domains  with  superlattice 
dimensions  (planes  ±1/3,  ±1.  'A  and  ±2/3,  ±2, ! A ).  Vacancy  rich  colums  are  frequently  seen  along 
[131]  and  [261],  Further  refined  simulations  using  supercell  to  describe  the  periodicity  of 
vacancies  is  required. 
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